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1. Summary 


A project to study the “Compressibility Effects on Dynamic Stall of Airfoils Undergo- 
ing Rapid Transient Pitching Motion” wets funded by AFOSR in FT -90. The objective of 
the study was to obtain a better understanding of unsteady flow separation as encountered 
in dynamic stall of airfoils, under local compressibility conditions and eventually control 
it for useful gains on real flight systems. 

The research was carried out in the Compressible Dynamic Stall Facility, CDSF , at 
the Fluid Mechanics Laboratory(FML) of NASA Ames Research Center. The facility 
can produce realistic nondimensional pitch rates experienced by fighter aircraft, which on 
model scale could be as high as 3600° /sec. Nonintrusive optical techniques were used 
for the measurements. The highlight of the effort was the development of a new real 
time interferometry method known as Point Diffraction Interferometry - PDI , for use in 
unsteady separated flows. This can yield instantaneous flow field density information (and 
hence pressure distributions in isentropic flows) over the airfoil. A key finding is that the 
dynamic stall vortex forms just as the airfoil leading edge separation bubble opens-up. A 
major result is the observation and quantification of multiple shocks over the airfoil near 
the leading edge. A quantitative analysis of the PDI images shows that pitching airfoils 
produce larger suction peaks than steady airfoils at the same Mach number prior to stall. 
The peak suction level reached just before stall develops is the same at all unsteady rates 
and decreases with increase in Mach number. The suction is lost once the dynamic stall 
vortex or vortical structure begins to convect. 

Based on the knowledge gained from this preliminary analysis of the data, efforts to 
control dynamic stall were initiated. The focus of this work was to axrive at a dynamically 
changing leading edge shape that produces only ‘acceptable’ airfoil pressure distributions 
over a large angle of attack range. 


2. Introduction 


Dynamic stall is an important problem and of particular interest to aircraft execut- 
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ing rapid maneuvers. The production of increased lift and stall delay during a transient 
maneuver could be effectively used to enhance the agility of fighter aircraft, provided the 
effect could be sustained for the duration of the maneuver. This requires a careful control 
of the process which is extremely complex 1,2 and depends on several factors such as Mach 
number, Reynolds number, pitch rate, airfoil geometry - leading edge shape - in particu- 
lar, transition etc. Among these, the compressibility effects (Mach number), leading edge 
boundary layer transition and leading edge shape (pressure distribution) play a more crit- 
ical role in altering the process. Compressibility effects set in at a very low free stream 
Mach number - 0.2 -0.3(McCroskey 1 , Chandrasekhara et al 3 , Chandrasekhara and Carr 4 ), 
and can cause premature stall. Production of shocks 3 over the airfoil and their interaction 
with the local boundary layer are partly responsible for this. However, the development 
of the leading edge adverse pressure gradient to the critical value at which flow separation 
occurs, at a lower angle of attack, is the main reason for the premature stall observed 
under compressibility conditions. 

The research effort described was aimed at understanding the underlying physical 
processes during the occurrence of dynamic stall of a 3in. chord NACA 0012 airfoil pitching 
transiently at rates of up to 3600° j sec in a free stream flow ranging from M = 0.2 - 0.45. 
A new nonintrusive optical flow measurement technique referred to as Point Diffraction 
Interferometry was developed with the funding to obtain real time quantitative flow field 
information. The dynamic motion of the airfoil at rapid rates precluded use of other 
conventional measurement methods. The report describes the results of the effort. 

3. Experimental Facility and Technique 


3.1. The Facility 

The experiments were performed in the Compressible Dynamic Stall Facility at NASA 
Ames Research Center. The CDSF 5 was established for conducting dynamic stall research 
in the FML. This facility is specifically designed for study of dynamic stall over a range of 
Mach numbers, using nonintrusive optical flow diagnostic techniques. It is operated as a 
part of the indraft tunnel complex at the FML. The CDSF is unique in that the airfoil is 


6 



supported between two lin. thick optical quality glass windows by pins that are smaller 
than the local airfoil thickness. Thus, the entire flow field including the airfoil surfaces can 
be viewed unobstructed by any support mechanism. This enables the study of the flow 
at the surface near the leading edge, where the dynamic stall vortex forms, as well as the 
flow field away from the airfoil. 

The rapid ramp- type pitching motion of the airfoil was produced by a hydraulic drive 
located on top of the test section. It can pitch an airfoil at constant pitch rates and is 
shown in Figure 1. 

The specifications of the hydraulic drive are as follows: 


angle of attack, a: 
pitch rate, d: 

maximum acceleration rate: 
change in a during acceleration: 
minimum acceleration time: 
free stream Mach number: 
airfoil chord: 

Reynolds number: 


0-60° 

0-3600 °/sec 
600,000 °/sec 2 
<6° of pitch 


4 ms 
0. 1-0.5 

3.0 in. 

180.000 - 840,000 


The pitch rate of 3600 °/sec on a 3in. chord airfoil corresponds to a 90°/sec pitching 
of a 10 ft. chord airplane wing at any given Mach number; thus, the rates obtainable 
from the design are directly applicable to flight conditions. In order to limit or isolate the 
effects of transients on separation, the change in angle of attack during acceleration and 
the acceleration time itself were limited to less than 6° and 4 ms, respectively. To properly 
simulate a maneuver, an angle of attack range of 0 - 60° was selected. To provide for 
reasonable experiment times, the facility has a recycle time of 2 seconds (30 runs/minute). 
The system uses both the airfoil position and velocity information in its feed back loops 
to properly perform any required maneuver which can be selected through software. The 
complete details of the final design are presented in Chandrasekhara and Carr 6 . 

The airfoil position was read by a digital optical encoder, whose output was input to 
the digital I/O board of a micro VAX II Work Station and timed with its internal clock. 





Figure 2 shows an example of the actual rates obtained, including the variation of the 
angle of attack during the various parts of the pitch-up motion in an experiment. It is 
clear that the airfoil pitches through its static stall angle at a constant rate of change of 
angle of attack. For the highest rate, the motion is completed in 15 ms, beyond which the 
system can be seen to be settling down (at the highest angle). All the tests were limited 
to the linear range. 

Since its original design, changes were made to the drive system to incorporate ad- 
ditional safety features and augment the protective features of the drive. These include 
installation of pressure gages to show the pressure in the system, pressure relief valves to 
bypass the flow directly to the tank - a precaution needed to mitigate the strong water 
hammer effects felt when the system was suddenly depressurized, adding new plumbing to 
control the pressure in the hydraulic accumulator during each pitch-up. electrical system 
modifications that de-energize the solenoid valve and c dump' the hydraulic line pressure, 
etc. Since the response of the feed back system would be altered by these components, a 
new system calibration was performed after tuning the feed back system properly to ensure 
stable operation of the system through the envelope of the experimental conditions. 

3.2. Instrumentation and Techniques 

The CDSF is equipped with a wide range of nonintrusive optical flow diagnostic instru- 
mentation such as stroboscopic schlieren, laser Doppler velocimetry, and point diffraction 
interferometry systems. Only the schlieren and the interferometry systems are discussed 
in this report. 

A. Stroboscopic schlieren studies 

A standard mirror based schlieren system was set up in a ; Z’ type configuration with 
a xenon arc lamp light source located at the focal length(10 ft.) of one of the mirrors. 
The beam passing through the test section was focused on to a vertical knife edge and 
imaging optics. The light source was triggered externally at the desired phase angles by 
an electronic circuit which compared the chosen phase angle of oscillation and the encoder 
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data from the drive system and produced a TTL pulse when a match occurred. No delays 
were found to be present between the events of matching the phase angle and the light 
flashing. 

B. Interferometry studies 

As stated earlier, the quantitative flow field density information was obtained using 
PDI. The PDI system used a CW/pulsed pumped Nd-YAG laser light source, with its 
beam expanded to 6in. to fill the entire field of view in the standard Z-type schlieren 
configuration. The optics were aligned to minimize astigmatism. A predeveloped, partially 
transmitting photographic plate was placed at the location of the knife edge as the PDI 
plate. Imaging optics were set up further downstream along the beam path for recording 
the flow. In operation, a pin hole was created in situ in the photographic plate with no flow 
in the test section. This acted as a point diffractor for the reference beam. Light deflected 
by the flow density changes (signal beam) focused to a slightly different spot, passed 
through the partially transmitting photographic plate and interfered with light passing 
through the pin-hole (which was the reference beam) to produce real time interference 
images, which were captured on Polaroid film. Fig. 3 shows a schematic of the final 
system used. 

Several different PDI plates were tried in an effort to improve the quality of the 
interferograms obtained. Some of these were: a green line filter, an unexposed holographic 
film, a predeveloped holographic plate etc. The trials showed that the best performance was 
obtained from a plate with an optical density of about 1.0. However, as is well known, film 
development is a function of several parameters and thus, trial and error is still involved. 
A measurement of the plate optical density is a good indication of the performance that 
can be obtained. 

One of the key parameters affecting the results of the technique is the size of the 
spot used as the point diffractor. A study to establish the sensitivity of the image quality 
and accuracy with respect to the spot size showed that spots of 50-80 \im diameter gave 
acceptable results over the whole flow field, including the dense region (with high frequency 
fringes) near the leading edge. These and other details have been documented in Ref. 7. 
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Ref. 7 and S describe more details of the technique and its implementation in the CDSF. 
C. Fringe Analysis and Image Processing 

A review of the available options for fringe analysis and image processing showed 
that the standard packages are very limited and almost exclusively, application specific. 
Thus, the task of developing a software package was undertaken in stages. In the present 
version, the package runs in a semi-automatic mode. It reads a digitized (256 gray level) 
interferogram and overlays the airfoil on the image. In addition, two slightly larger airfoils 
are also drawn. The user interactively picks the intersection of the fringes with the airfoil 
upper and lower surfaces separately. When the fringe density is high, it permits the user to 
go into the ‘off-body’ mode and pick fringes along a line parallel to and away from the airfoil 
surface (on one of the larger airfoils) where the fringes are farther apart. These are then 
transferred to the airfoil surface. Typically, the fringes near the leading edge region are very 
dense in the flow field being studied owing to the very large local density gradients. Optical 
noise introduced by the shadowgraph effect in this region lowers the contrast, making it 
an ideal location where the off-body mode needs to be invoked. The package provides an 
output data file containing the various physical variables, using isentropic flow relations in 
a format suitable for the FML/NASA ARC standard QPLOT graphics package. Presently, 
the data is being processed on an IRIS 4000 Series Work Station and it takes about 3-5 
minutes per image for processing. 

4. Results and Discussion 

The results are presented in the form of schlieren pictures, PDI images and plots of 
data retrieved from these flow images. The schlieren results are described in ref. 3 and are 
simply reproduced here for completeness. 

4.1. Effect of Mach Number 

Fig. 4 compares the schlieren pictures at different Mach numbers for a + = 0.03 and 
a = 17°. It can be seen that for M < 0.3, the vortex is at about 50% chord location. 
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In addition, the vertical extent of the flow is nearly the same for M = 0.2, 0.25 and 0.3. 
However, for M > 0.3. the dynamic stall vortex moves successively closer to the trailing 
edge and the flow scales can be seen to have increased as well. Movement of the vortex 
downstream indicates that the flow is approaching the deep stall state and thus, it is 
clear from the figure that as the Mach number is increased above 0.3, deep stall occurs at 
progressively lower angles of attack. 

Fig. 5a shows the effect of Mach number on dynamic stall for a pitch rate, o + = 
0.025, and the corresponding results for = 0.035 are shown in Fig. 5b. Plotted in it 
are the successive locations of the center of the dynamic stall vortex as a function of the 
instantaneous angle of attack at different Mach numbers. It can be seen in both figures 
that the vortex appears at lower angles of attack as the Mach number increases. This also 
leads to the result that the vortex moves past the trailing edge at lower angles of attack for 
higher Mach numbers, causing deep dynamic stall to occur earlier in the pitch-up motion. 
Significant decrease in the angle of attack occurs for the same x/c location for M > 0.3 
and thus, M = 0.3 can be considered to be the limit when compressibility effects set in. 
Consider for example Fig. 5a, for x/c as 0.6, the center of the vortex is at a = 16.5° for 
M = 0.3, and a = 14° for M = 0.45. Similarly, in Fig. 5b, the vortex is at 60% chord 
location at a = 19° for M = 0.3; at M = 0.4, the corresponding angle of attack = 17.2°. 
Similar results were obtained at other pitch rates. 

4.2. Effect of Pitch Rate 

Fig. 6a through 6d show the vortex center locations over the airfoil plotted as a 
function of the angle of attack at different pitch rates for M = 0.2, 0.35, 0.4 and 0.45 
respectively. It can be seen in all the figures that the vortex is retained on the surface of 
the airfoil to higher angles of attack as the pitch rate is increased. The trend is monotonic 
with increasing pitch rate. For example at M = 0.45, the vortex is on the surface even at 
a = 18° at = 0.03, whereas the static stall angle for this case is ^ 9.5° as determined 
from the schlieren images. For = 0.020, deep dynamic stall occurs at a = 15.5°. For 
M = 0.35, the deep stall angle is « 23° for a+ = 0.04, and the static stall angle is 11.6°. 
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4.3. Leading Edge Supersonic Flow 

There has been discussion for a long time about the possibility of a shock developing 
in the leading edge flow. Fig. 7 provides conclusive proof that a shock, in fact, multiple 
shocks form in the flow. It shows a schlieren picture of the flow at M = 0.45, a — 12.6° 
at a + = 0.0313 and a PDI image for the same conditions. The latter actually allows 
quantification of the flow. A simple fringe counting shows that the local Mach number 
ahead of the shock is greater than 1.0 and at the foot of the shock, it is about 1.2. Although 
the flow is only weakly supersonic, the shock causes the leading edge laminar boundary 
layer to separate. This separated free shear layer develops a waviness, which causes the 
flow downstream of the shock to go through a series of accelerations and decelerations. As 
the flow negotiates the crests and valleys of this wavy shear layer, expansion waves and 
compression waves develop, causing a series of shocks. The last shock in the series appears 
to be the strongest and the flow becomes subsonic downstream. Further analysis of this 
interesting basic flow is ongoing. The shocks are repeatable, and have been found to be 
present over an angle of attack range of about one degree. 

4.4. Separation Bubble and Dynamic Stall 

Fig. 8a is a point diffraction interferogram of the flow at M = 0.3, a = 12° and 
a + = 0.03. This image reveals some important features of the flow. Firstly, the dark 
fringe on the lower surface slightly aft of the leading edge is the stagnation point fringe. 
The suction pressure developed by the airfoil causes the local flow to accelerate, resulting 
in strong density changes, which is seen in the figure as a concentration of high frequency 
fringes. The close spacing of the fringes also means that the flow gradients are very high. 
In fact, 21 dark fringes are present in 1 millimeter in this image, indicating that the local 
Mach number is 0.71 and the local pressure coefficient is -3.75 at x/c = 0.01. Downstream 
of this the steep adverse pressure gradient region where the classic laminar flow separation 
occurs. The separated shear layer reattaches after it transitions into a turbulent layer, 
forming a bubble. The accompanying pressure distribution obtained by processing the 
image using the fringe analysis software developed is plotted in Fig. 8b. It shows the 
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suction peak, the drop in suction due to the adverse pressure following it and the laminar 
separation bubble, the later is indicated by the plateau in the distribution. Downstream 
of the plateau, the pressure gradually recovers to the free stream value. When several such 
pictures were analyzed, it was found that the nondimensional pressure gradient ^ at 

which the bubble opens up is about 100. As the shear layer moves outwards, the dynamic 
stall vortex forms. During this stage, the airfoil continues to develop additional suction 
(since the shear layer is still attached at the leading edge). However, once the vortex grows 
and begins to convect, the suction drops. Interestingly, during the convection process, the 
pressure distribution flattens out, leading to continued production of lift. These details 
are discussed further in ref. 9. Only sample pressures distributions are shown in Fig. 9 
for angles of attack ranging from 11-15 degrees. 

4.5. Leading Edge Flow Development vs. Unsteadiness 

The maximum suction pressure attained in the flow just before the vortex formed and 
became discernible in the PDI pictures was estimated by fringe counting. Its variation 
with the pitch rate is shown in Fig. 10 9 . Included in it are the steady flow values when 
stall occurred. A significant result is the near constancy of the Cp m - n at a given Mach 
number for all unsteady flow cases plotted. The values are always higher than the steady 
flow levels, clearly indicating the fact that transient pitching of the airfoil keeps the fluid 
attached enabling the airfoil to develop larger suction. Further, since the airfoil reaches a 
higher angle of attack before the vortex forms, stall delay is concomitant with the pitching 
process. It can be inferred from this that the leading edge pressure distribution is modified 
considerably by pitching the airfoil. The striking result in the figure is the strong depen- 
dence of the peak suction pressure coefficient on Mach number, decreasing monotonically 
with increasing Mach number. 

Fig. 11 compares the development of the peak suction pressure coefficient with angle 
of attack for M = 0.4 at different degrees of unsteadiness. It is evident that the pressure 
development lags that of steady flow by a few degrees, depending on the rate of pitching. 
In addition, the unsteady flow also can reach a higher suction peak of -4.5 (and hence 
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higher lift) as opposed to the steady flow value of about -3.0. The peak remains high as 
the vortex forms and grows, but once vortex convection starts, it begins to drop. The 
vortex formation and growth occurs over an angle of attack range of about half-a-degree. 

4.6. Control of Dynamic Stall 

The PDI studies have revealed that the leading edge pressure distribution develop- 
ment is delayed in unsteady flows. To use dynamic stall benefits in aircraft systems, this 
delay needs to be sustained and even, increased further. Since the dynamic stall vortex 
occurrence correlates with the C Prnin values over the leading edge (presumably the level 
at which the critical adverse pressure gradient for stall is encountered downstream) a pos- 
sibility exists that if the pressure distribution is altered suitably, stall control could be 
achieved. Towards this goal, a flexible leading edge over an airfoil offers excellent poten- 
tial, since the airfoil geometry could be instantaneously modified to produce a favorable 
pressure distribution always . A material developed and patented at NASA ARC has been 
found to have the necessary properties for this attempt. This material is a sheet of plastic, 
in which are embedded two strips of a conducting material. When subjected to a like 
charge, the repulsive force between the two strips produces sizeable surface deformation. 
Such a device is now being used to de-ice aircraft wings. 

In the current application, the idea is to replace the airfoil leading edge with this 
material, apply a controlled voltage to it and produce a dynamically deformable leading 
edge (DDLE). An existing airfoil made from this sheet was subjected to different voltages. 
The rapid surface deformation cycle that resulted lasted about 3 milliseconds. Images of 
the surface were recorded with a high speed IMA-CON camera (at rates of up to 25,000 
frames/sec). The camera imaged up to 16 frames at the selected rate. The challenge was 
to identify the portion of the surface that was being studied and illuminate it properly for 
good photographs to be obtained. An orthogonal grid of retro-reflecting tape glued to the 
airfoil and lit by strobe lights provided this information. The camera and the strobe lights 
were triggered by the same pulse generator driving the leading edge deformation circuit. 
The camera/ strobe circuit incorporated a controllable delay so that any desired portion of 
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the deformation cycle could be recorded. Fig. 12 shows the schematic of the arrangement 
used. The resulting images were digitized and processed on an IRIS Work Station, using 
software developed in-house. The deflections were measured from a stationary reference. 
A sample plot of the deformation is shown in Fig. 13, where a section of the surface is 
plotted at three instants of time. 

The deformation of the strips used was found to be three dimensional which is obvi- 
ously unacceptable for dynamic stall control. Design effort is now underway to achieve a 
two-dimensional surface deflection through a deformation cycle. The aerodynamic behav- 
ior of the airfoil will then be examined to determine its effectiveness in overcoming both 
static and dynamic stall. 


5. Concluding Remarks 

A comprehensive study of dynamic stall of a transiently pitching airfoil has been 
conducted. It has shown that: 

1. Compressibility effects set in at a Mach number of 0.3. 

2. Dynamic stall arises out of the bursting of a laminar separation bubble near the 
leading edge on the upper surface of the airfoil. 

3. Multiple shocks form over the airfoil at M = 0.4 and above, for certain angles of 
attack. The associated fluid flow interactions are very complicated. 

4. The suction pressure coefficient is a function of the Mach number only and decreases 
with increasing Mach number. 

5. The delay of stall is brought about by a delay in the development of the leading edge 
pressure distribution, when compared to steady flow under the same conditions. Increasing 
unsteadiness increases the delay. 

6. The correlation between the airfoil peak suction pressure coefficient and dynamic 
stall occurrence leads to stall control ideas using adaptive geometries, which are being 
attempted. The challenge is to form a flexible leading edge airfoil that can deflect only in 
two dimensions. The effects of changing leading edge geometry on fluid flow still needs to 
be investigated. 
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2.5 GALLON 
/ACCUMULATOR 



SERVO VALVE 1 



Fig. 1. Diagram of the In-draft Wind Tunnel with the Constant Pitch Rate 
Mechanism 


U = 0.45, a~ « 0.03 , SLOPE 0-57 DEC = 3507.0 

U - 0.35. a + = 0.025, SLOPE 0-57 DEG = 2256.0 

M - 0.25. a~ - 0.02 . SLOPE 0-57 DEC - 1233.0 



time, sec 


Fig. 2. Time-history of Ramp for Sample Test Conditions 
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f- BEAM FROM 
LASER 



Fig. 3. Schematic of the Point Diffraction Interferometer. 



M = 0.35 M = 0.40 M = 0.45 


Fig. 4. Effect of Mach Number on Dynamic Stall of a Rapidly Pitching Airfoil: o + = 
0.03, a = 17deg. 
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0.025, (b) a+ = 0.035. 







tative Effects 





(a) Schlieren Picture 



(b) Point Diffraction Interferogram 



Fig. 7. Multiple Shocks on a Rapidly Pitching Airfoil: M = 0.45, a + = 0.0313, a — 12.6°. 


22 








Fig. 9. Pressure Distributions Over a Transiently Pitching Airfoil During Dynamic Stall 
Onset: M = 0.3, a + = 0.04. 
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Pitching Airfoil. 
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